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I.   INTRODUCTION 

The  field  of  Signal  Processing  is  expanding  rapidly, 
indicating  an  increasing  importance.  Today,  digital  signal 
processing  is  being  extensively  explored  and  employed. 
High  speed  digital  computers  with  necessary  peripheral 
equipment  such  as  analog-to-digital  and  digital-to-analog 
converters  are  required.  Traditional  analog  systems  as  radar 
and  sonar  are  going  digital.  There  exists  a  wealth  of  other 
digital  signal  processing  applications. 

For  some  applications,  however,  the  use  of  digital 
processing  may  be  termed  "overkill",  even  if  there  is  a 
distinct  requirement  for  the  inherent  advantages  of  digital 
processing  over  analog  processing. 

A  fairly  new  device  technology,  the  Charge  Transport 
Device  technology  (CTD's)  offers  new  ways  of  signal  processing 
which  could  be  the  answer  for  such  cases. 

Charge  Transport  Devices  refer  to  several  families  of  devices: 
Bucket  Brigade  Devices  (BBD's),  Charge  Coupled  Devices  (CCD's), 
and  Reticon  Devices  (i.e.  SAD-100  and  TAD-12)  which  are 
described  in  reference  {3}. 

The  unique  feature  of  CTD's  is  that  the  signal  is 
represented  by  discrete  charge  packets,  which  are  moved 
through  the  device  at  a  speed  determined  by  the  clock 
frequency.  The  delay  of  the  device  is  thus  controlled  by  the 
clock  frequency,  and  may  be  easily  and  momentarily  changed. 
An  excellent  review  of  how  the  CCD  works  is  given  in 
reference  { 8 } . 

CTD's  find  wide  applications  in  imaging  (both  low  light, 
normal  light  and  infrared),  analog  time  delay,  multiplexing, 
transversal  filtering,  with  fixed  and  variable  weighting, 
recursive  filtering,  correlators,  chirp-z  transformers  and 
so  on.     The  general  advantages  of  CTD  signal  processors 
over  conventional  digital  processing  are  lower  cost,  power, 
weight  and  size. 
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Comb  filters  are  described  by  their  periodic  transfer 
characteristics  in  the  frequency  domain.  They  can  be  classi- 
fied into  two  general  types.  The  first  type  is  the  bandstcp 
or  canceller  type  shown  in  Table  II.   It  has  strong  attenu- 
ation in  a  narrow  neighbourhood  of  a  series  of  periodically 
separated  frequencies  and  good  transmission  in  between. 
The  second  type  is  the  bandpass,  or  integrator  type  (Table  II) 
It  has  good  transmission  in  a  narrow  neighbourhood  of  a 
series  of  periodically  separated  frequencies  and  strong 
attenuation  in  between.   All  possible  implementation 
techniques  for  such  comb  filters  use  some  kind  of  delay  de- 
vice . 

In  the  analog  case,  a  quartz  delay  line  has  been 
extensively  used.  It  is,  however,  temperature  dependent 
and  results  have  in  general  not  been  very  satisfactory. 
For  the  digital  case  the  delay  device  is  the  shift  register, 
in  which  the  amount  of  delay  is  determined  by  the  employed 
clock  frequency. 

Today  a  new  family  of  comb  filters  is  being  developed,  using 
the  new  CTD's.  In  these,  analog  signals  are  first  sampled 
and  then  delayed.  As  a  result  of  the  sampling,  these  devices 
not  only  have  standard  analog  properties  but  also  some 
digital  properties  such  as  aliasing  and  stability  of  delay. 
A  CTD  consists  of  N  serial  delay  stages.  A  signal  will  then 
be  delayed  N  clock  periods  before  it  is  fed  back  to  the 
input  in  the  recursive  case.  The  frequency  of  recursive 
operation,  f  ,  is  then  only   l;f  ,  where  f   is  the  clock 
frequency,  or  sampling  frequency.  Therefore,  in  radar  terms, 
f   is  equal  to  the  Pulse  Repetition  Frequency  (PRF) ,  and 

there  are  N  teeth  of  the  comb  within  the  Nyquist  frequency 

2 
range  from  0  to  f  /2. 

So  far  recursive  filtering  using  CCD's  have  not  been 
extensively  studied  in  the  industry,  since  it  was  believed 
that  it  can  only  be  done  with  extensive  external  electronic 

support  such  as  fast  and  accurate  operational  amplifiers  etc. 

{5}. 
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In  fact,  reference  {4}  was  the  only  paper  on  CCD-  recursive 

filters   presented  at  the  recent  conference  on  CCD-appli- 

cations  in  San  Diego. 

This  thesis  will  study  the  second  order  recursive  comb  filter 

and  employ  a  simulation  scheme  to  study  its  application 

as  a  canceller  in  an  MTI-radar. 
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II.   SAMPLED  ANALOG  SECOND  ORDER  RECURSIVE  COME  FILTERS 

A.    THEORY 

1 .   Analog  Filter 

Since  the  sampled  analog  recursive  filters  are 
implemented  by  the  same  circuit  configuration  used  for  digital 
recursive  filters,  we  may  apply  the  existing  well-developed 
theory  for  such  filters.  However,  proper  modifications  must 
be  made  to  account  for  the  fact  that  the  signal  is  now  sampled 
analog  and  that  there  now  are  N  delay  stages  instead  of  one. 
It  is  known  that  two  general  design  approaches  have  been 
developed  for  the  IIR  (infinite  impulse  response)  type  digital 
recursive  filter,  namely  the  indirect  and  direct  methods. 
The  indirect  method  starts  with  the  transfer  function  in  the 
LaPlace  transform  variable   s   (H(s)),  of  a  continuous 
filter  and  follows  by  performing  another  transform  to  obtain 
the  corresponding  transfer  function  in  the  discrete  time 
variable   z  .   Several  different  transforms  are  available 
for  this  purpose,  for  example  mapping  by  differential 
transform,  impulse  invariant  (standard  z  )  transform, 
bilinear  z- transform  and  matched  z- transform. 
The  direct  method  is  digital  design  in  frequency  or  time 
domain  using  some  type  of  computer  aided  design  procedure. 

Here  the  indirect  approach  will  be  used: 
Given  a  prescribed  response,  there  are  several  methods  by 
which  one  can  develop  a  transfer  function  in  the  variable  s. 
Examples  are:   Butterworth,  Chebyshev  and  elliptical  filters 
etc{l}.  In  this  thesis  the  Butterworth  approach  is  used  to 
obtain  the  transfer  function  in  the  s-domain. 

a.   For  a  Low  Pass  Filter 

For  the  Butterworth  filter  in  the  second  order 
case  the  equation  is  given  as: 

H(ju))2  r 


1  +  W 
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where  cu   is  the  3-db  cutoff  frequency. 
Converting  to  s-domain  and  expanding: 

H(s)-  H(-s)  =         1 


2N  Let  N  =2 

s 

J  x 

By  evaluating  the  poles  in  the  s-plane  and  discarding  the 
two  poles  in  the  right  hand  plane,  one  gets: 

H(s)  =         1 


2  2 

s   +  2u  s  +  u 

X        X 


This  is  the  desired  transfer  function  in  the  s-domain. 
b.   For  High  Pass  Filter 

A  similar  derivation  for  the  high  pass  case 
gives : 

H(s)  =  s2 


2  2 

s   +  2w  s   +   a) 

X  X 


2  .   Transformation  to  Z-dornain 

a.  Low  Pass  Filter 

By  indirect  design  method,  using  the  bilinear 
z-transf orm: 

s  =  2(z-l) 
TTz+lT 

one   gets  the  following  transfer  function  in  z: 

2 

CO 

H(z)    =  X 


\T    z  +  l] 


A      2    z-1    0  2 

T    z+1        x  x 


By  employing  a  little  algebra  one  will  end  up  with  this 
equation: 
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H(z)    = 


K    (    aQ    +    a1    z        +    a2    z~2) 


1   +   b1    z"1    +      b2    z"2 


where:                                 „2    ,     2 
K    =      -*- 


4    +    2/2    0)      T    +    w    2    T2 

X  X 


a      =    1 
o 


ai  =  2 


a      =    1 

2         2 
2    a)  T 

k      -             x 
bl    "      " —— 


4    +    2/2       03      T+      co2T2 

X  X 

4    +      w    2    T2    -    2V7!  w   T 

K         -  X  X 

b2    ~      ' — " — 

4    +    2    /2oj   T      +      u)    2    T2 
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b.   High  Pass  Filter 

Starting  from  the  s-domain  equation 


H(s)  =   -T— -S y 

S   +  /2  03   s   +0) 
X  X 


and  using  the  bilinear  transform  again,  one  end  up  with 

K(   a      +   a.    z        +   a„   z      ' 

H(z)    =      2- L_ 2 


-1  -2 

1    +   b,    z        +    b^    z 


where : 


K   = 


H    +    2    /7  u   T   +      co    2    T2 

X  X 


a      =    an      =    1 
o  2 


ai  =  -2 


bi 


2    2 
2    co      T 
x 


b2 


4    +    2    /2    03   T      +      co    2T2 

X  X 


4      -2/203T  +      03      T 
x  x 

4    +       2    /2    03    2T      +      co    2T2 

X  X 


It  is  noted  that  oj  and  T  always  occur  together  in  the  final 
expressions  and  that  there  are  only  two  differences  between 
the  high  pass  and  low  pass  cases:  The  size  of  K,  and  the  sign 
of  a,.   But  K  is  merely  a  scaling  factor,  and  has  no  influence 
upon  the  shape  of  the  frequency  response. 
This  means  that  when  the  coefficients  are  kept  constant,  the 
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TABLE 


GO     T 
X 

bl 

u>   T 

X 

b2 

0 

-2. 

0 

1. 

.1 

-1.86 

.5 

.5 

.5 

-1.32 

1. 

.28 

.7 

-1.09 

1.5 

.18 

1. 

-.77 

2. 

.17 

1.3 

-.49 

3. 

.21 

1.6 

-.26 

4. 

.28 

1.9 

-.06 

5. 

.34 

2. 

0 

6. 

.4 

3. 

.47 

7. 

.46 

4. 

.77 

6. 

1.24 

Region  I 


Region     II 


COJ 
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cutoff  frequency  will  change  by  varying  T. 

Table  I  and  Figure  2.1  shows  how  b,  and  b0  varies  with  to  T. 

Table  II  is  a  summary  of  how  the  frequency  response  behave 
as  w^T  is  varied,  giving  the  ranges  where  the  filter  acts  as 
integrator  or  canceller. 

Since  region  I  (Figure  2.1)  is  the  area  of  interest  in  this 
thesis,  the  experimental  study  was  carried  out  here. 

3.   Direct  Design  Method  in  Z-Domain 

a.   By  Trial  and  Error  Method 

By  carrying  out  the  design  in  the  z-domain, 
transformations  from  HCs)  to  H(z)  can  be  avoided: 


Z-plane 


FIGURE   2.2.   Z-Domain  Unit  Circle 

Figure  2.2  shows  the  z-plane,  where  the  entire  left  half  of 
the  s-plane  has  been  transformed  into  the  inside  of  the  unit 
circle.   By  moving  the  two  poles  and  two  zeros  (second  order) 
around  and  using  the  geometric  method  {1}  to  calculate  the 
amplitude-  frequency  response,  it  is  possible  to  find  a 
combination  of  the  singularities  that  gives  a  response  quite 
close  to  the-  desired  response.   Appendix  A  shows  the  result 
of  such  a  study. 

b.   By  Computer  Approach 

Several  computer  programs  {1}  have  been 
developed  that  will  optimize  the  search  of  coefficients, 
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TABLE   II 


Relationship  between  w _T  and  Filtertype  for  Highpass 
and  lowpass  filters  using  the  bilinear  z-transfonrt  method 


Filter  Type 


OJ     T 
x 


Filter  Shape 


High  Pass 


Low  Pass 


<2 

(Region  I, Canceller) 


>2 

(Region  II , Integrator ) 


<2 


(Region  I,  Integrator) 


>2 


(Region  II,  Canceller) 


s  =  sampling  frequency     N=  number  of  delay  stages 
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based  on  different  approaches  in  reducing  the  deviation 
between  presented  response  and  filter  response. 
This  approach  was  not  employed  here. 

B.    EXPERIMENT 

1.    System  Configuration 


FIGURE  2.3.   Second  Order  Recursive  Filter 

Figure  2.3  shows  the  block  diagram  of  the  canonical  form  of 
a  second  order  recursive  filter.  Serial  analog  delay  (SAD) 
from  Reticcn  Co  were  used  as  delay  lines.   20  K  potentio- 
meters were  used  to  implement  the  coefficients.  In  the  cases 
where  coefficients  greater  than  one  were  needed  for  the  "a"- 
coef f icients ,  the  numerator  in  the  transfer  function  H(z)  was 
scaled  down  so  that  the  "a"-coeff icients  were  normalized. 
Several  methods  to  implement  the  coefficients  using  the 
potentiometers  are  discussed  in  references  {2}  and  {3}. 
In  this  work  Figure  2.M  is  presented,  giving  the  calibration 
of  potentiometer  setting  for  desired  coefficient  value.  This 
curve  was  obtained  by  simply  measuring  input  and  output 
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across  the  potentiometer. 

2.   Results 

The  clock  frequency  for  the  SAD-100  (Reticon  Co) 
device  was  72  KHZ.  Since  the  effective  clock  frequency  was 
four  times  72  KHZ,  and  there  was  96  stages  within  the  device, 
the  corresponding  nulls  of  a  high  pass  comb  filter  should 
occur  at  multiples  of: 


7  2  KHZ 


96 


3  KHZ 


A  computer  program  (APPENDIX  B)  was  used  to  calculate  and 
plot  the  theoretical  frequency  response  curve, and  also  to 
superimpose  the  actual  data  points  on  the  same  plot. 
The  "a"-  coefficients  were  kept  fixed,  and  several  combina- 
tions of  b-,  and  b,-,  were  used.   They  were  selected  from 
Figure  2.1,  and  b,  was  negative  in  all  cases  studied  in  this 
thesis . 

By  using  the  values  of  b-,  and  b~  corresponding  to  to  T ,  the 
resulting  response  should  be  flat-topped  in  every  case. 
The  implemented  coefficients  were: 


High  Pass 
(Canceller) 


Low  Pass 
(Integrator) 


a 
o 

al 

a2 

1 

bl 

b2 

1 

.5 

-1.0 

.  5 

-1.0 

.38 

r 

.  a 

-1.0 

.  5 

-0.7 

.27 

.5 

-1.0 

.5 

-0.4 

.21 

.  5 

-1.0 

.5 

-0.3 

.2 

.5 

-1.0 

.  5 

-0.2 

.19 

.5 

1.0 

.5 

-1.0 

.38 

.5 

1.0 

.5 

-0.7 

.27 

.5 

1.0 

.5 

-0.4 

.21 

)          .5 

1.0 

.5 

-0.3 

.2 

.5 

1.0 

.5 

i — , — 

-0.2 

.19 
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The  coefficients  b1  and  b^   were  all  picked  in  Region  I 
(Figure  2.1). 

Figures  2.9  through  2.18  show  the  series  of  results. 
There  seemed  to  be  the  same  problem  with  the  second  order 
filter  as  was  experienced  with  the  first  order  filter{2}, 
namely  a  skewness  at  higher  frequencies.  Also  the  results 
show  an  undesired  dip  in  the  passbands ,  indicating  that  the 
coefficients  are  not  correctly  implemented. 
Suspecting  that  the  experimental  implementation  inherently 
makes  the  system  frequency  dependent,  several  checks  were 
made  for  its  investigation.  The  gain  of  individual  delay  lines 
was  checked  for  frequency  dependence,  using  both  the  oscillo- 
scope and  a  narrowband  voltmeter. 

The  individual  potentiometers  input-output  characteristics 
were  measured  when  the  filter  loop  was  closed. 
It  was  found  that  the  gains  of  the  delay  lines  were  slightly 
frequency  dependent  as  measured  on  the  oscilloscope,  and 
severely  and  inconsistently  frequency  dependent  when 
measured  on  the  narrowband  voltmeter,  indicating  loading 
effects  in  the  latter  case;  (the  input  impedance  of  the 
voltmeter  was  100  Kohm) . 

Figures  2.5  through  2.8  show  the  calibration  of  the 
potentiometer  settings  including  the  frequency  dependence 
for  the  coefficients  used  in  the  previous  measurements. 
Figures  2. IS  through  2.21  show  the  frequency  response  when 
the  frequency  compensated  coefficients  are  used. 
The  potentiometer  settings  were  changed  every  3  KHZ. 
During  several  of  the  "runs"  it  can  be  seen  that  there  is 
a  slight  shift  in  the  position  of  the  experimental  curve 
relative  to  the  theoretical  curve.  This  is  due  to 
uncontrollable  drift  of  the  basic  clock  frequency  in  the 
course  of  the  datataking. 

The  very  latest  set  of  data  (Figure  2.21),  was  taken  while 
adjusting  a   and  a2  at  every  null,  after  first  having  set  all 
coefficients  according  to  the  curves  given  in  Figures  2.5  - 


25 


2.8,  at  that  null. 

Table  III  shows  the  potentiometer • settings  for  a  =  0.5, 

a,=  -1.0,  a2=  0.5,  b1=  -0.7,  b2=  0.27: 


TABLE   III 


From  Curves 

After  readjusting 

Frequency 

a 
o 

al 

a„ 

z 

bl 

D2 

a 
o 

a  r\ 

3  KHZ 

6   " 

9   " 

12   " 

15   " 

18   " 

.714 

.714 
.714 
.715 
.718 
.722 

-1.0 
-1.0 
-1.0 
-1.0 
-1.0 
-1.0 

.846 

.846 

.847 

.848 

.85 

.852 

-.873 

-.  873 

-.875 

-.876 

-0.877 

-.877 

.407 

.405 

.403 

.4 

.398 

.  396 

.652 
.652 
.689 
.685 
.685 
.68 

.885 
.885 
.863 
.859 
.859 
.862 

Figures  2.19  through  2.21  do  not  show  any  sign  of  the 
skewness  of  previous  data  sets,  but  the  shape  of  the 
observed  curves  still  deviate  significantly  from   the 
theoretical  curves,  indicating  that  the  proper  coefficients 
have  not  been  found  for  all  frequencies. 


26 


i4 


IS1 
-■  cm  K 


.Q 


o 


00 
CM 


CN 


CN 


CN 


O 
CN 


00 


..   VO 


O 

rH 


CO 


KD 


"vT 


CO 


^9S 

^oa 

CO 
CO 

CO 

CO 

in 

CO 


CO 


in 


r- 


27 


o 
II 


-**- 


^8S      ^Od 

— — ^r -ft- 


^r 


CM 


II 
CM 

-a 
it 


ii 

H 

.Q 

4h 
O 

c 

O 
•H 
-H 

•H 

& 

rti 
> 

O 

G 

0) 

D^ 
0) 
& 


CD 
CN 


CD 
H 


CN 


co 

-'     CM 


£uyq.q.as  ■  ;oa 


<- 


V0 


CM 

rH 

yO 

«X) 

O 

o 

cn 
in 


CO 


LD 

O 


to 


CN 


O 
CN 


O 


CM 

H 


■■  CO 


<^ 


00 
CO 


X) 


H 
,Q 

<H 
O 

o 

•H 
•H 

rd 
•H 

> 

o 
r. 

0) 

O 


CN 

w 

CD 
M 


29 


o 

CO 


CO 
CM 


CM 


..     II 

CM 


o 

CM 


J.  CnJ 


CD 


II 

H 
4-rQ 


o 


CM 


oo  4- 


^"  .. 


CO 


en 

CM 


— +- 

OO 
CM 


CM 


CD 
CM 


CM 


CM 


CO 
CM 


CM 
CM 


CM 


O 
CM 


CD 

H 


30 


+*  t-*-* 


1 — ++^±±+1-.^ 


'\ 


*  ^tt:^ 


:+  +    +_+ 


+ 


_J-_ 


±=±=3=_ 


BO 


X^E^Z^n 


o         co         vo         -*  cv         o 

•         *         ♦         •  • 

H             O             O            O  O 

FIGURE  2.9a.  Theoretical  and  Observed  Frequency  Response  of  High  Pass 
Filter;  t^  =  -1.0  ,  t>2  =  0.38 
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Magnitude  (normalized) 

FIGURE  2,9b.  Theoretic*.!  and  Observed  Response  of  High  Pass  Filter; 
b  -  -1.0  ,  b2  -  0.3S 
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Magnitude  (nonralized) 


FIGURE  2.10a.  Theoretical  and  Observed  Response  of  Low  Pass  Filter; 
h  «  1.0  ,  b„  =  0.38 
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Magnitude  (normalized) 


FIGURE  2. 10b.  Theoretical  and  Observed  Response  of  Low  Pass  Filter: 
b  -  1.0,  b2  -  0.38 
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FIGURE  2.11a     Theoretical  and  Observed  Response  of  High  Pass  Kilter; 
bl  =  ~°'7'    b2  =  °'27 
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Magnitude  (normalized) 


FIGURE  2.  lib.  Theoretical  and  Observed  Response  of  High  Pass  Filter; 
b1  =  -0.7,  b2  =  0.2? 
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FIGURE  2.12b.  Theoretical  and  Observed  Response  of  Low  Pass  Filter j 
b1  =  °*7'   b2  =  °'27 
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FIGURE  2.14a.  Theoretical  and  Observed  Response  of  Low  Pass  Filter; 
b  -  0.4,  b2  =  0.21 
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FIGURE  2.14b.  Theoretical  and  Observed  Response  of  Low  Pass  Filter; 
b1  =  O.Ut     b2=  0.21 
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FIGURE  2.15b.  Theoretical  and  Observed  Response  of  High  Pass  Filter; 
b,  =  -  0.3,  b0  -  0.2 
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FIGURK  ?..17b.  Theoretical  and  Observed  Response  of  High  Pass  Filter; 
b1  -  -  0.2,  b2  -  0.19 
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FIGURE  2.18b.    Theoretical  and  Observed    Response  of  Low  Pass  Filter; 
b1  -  0.2,       b0  =  0.19 
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FIGURE  2.19a.   Theoretical  and  Observed  Response  of  High  pass  Filter; 
b,  -  -1.0;     b0  =  0.38 
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FIGURE  2.2  0a.  Theoretical  and  Observed  Re 
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FIGURE  2.21.    Theoretical  and  Observed  Response  of  His*  Pass  Filter; 
b    =  .0.7,     b9  -     0.27 
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III.   APPLICATION  TO  MTI-RADAR 

A.    THEORY 

1 .   Statement  of  the  Radar  Signal  Processing  Requirement 

The  majority  of  existing  radars  are  pulsed  radar 

systems.  Thus  the  radar  return  is  in  effect  a  sampled  signal, 

and  the  sampling  width  is  the  pulse  width  (PW)  of  the  radar. 

Figure  3.1a.  is  a  time-amplitude  plot  of  radar  signals. 

At  multiples  of  T=  t^v,  the  transmitted  pulses  are  shown. 

The  time  interval  between  each  transmitted  pulse,  called 

the  Pulse  Repetition  Interval  (PRI) ,  target  returns  will 

appear.   In  Figure  3.1a  two  stationary  targets  and  one  moving 

target  are  shown.   The  stationary  target  return  may  be  either 

positive  or  negative  depending  upon  the  phase  of  the  returned 

signal.  There  is,  however,  no  change  from  interval  to  interval, 

a  fact  which  is  exploited  in  the  MTI  canceller. 

The  moving  target  however,  is  changing  in  amplitude  every 

interval,  at  a  rate  which  is  the  doppler  frequency. 

Figure  3.1b  shows  how  the  pulse  train  of  Figure  3.1a 

would  look  as  presented  on  an  oscilloscope  or  A-scope. 

Here  every  PRI  is  shown  on  top  of  each  other,  and  the  moving 

target  is  exibiting  the  socalled  "Butterfly  effect"  {6} 

while  the  stationary  targets  remain  fixed  in  amplitude  and 

time. 

To  get  the  frequency  representation  of  a  periodic  pulse 

train  in  the  time  domain,  one  can  take  the  Fourier  transform. 

The  result  would  be  a    sin  x   -type  waveshape.   If  only  one 

x 
pulse  were  transmitted,  the  frequency  spectrum  would  be 

continuous,  and  the  first  null  of  the  curve  would  appear  at 

1 
"  PW*  But  when  transmitting  a  train  of  pulses  spaced  in  time 

at  T  -  ttttpt  ,  ideally  each  target  will  return  a  square  pulse 

at  PRI  equal  to  T.   The  frequency  spCctrum  of  this  periodic 
signal  can  be  found  from  the  Fourier  series  if  the  pulse 
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train  is  infinite: 


jnco  t 

f(t)  =  E  C   e    -1 

n 

—  00 


where 


1      1/2  -jno^t 

2 


—^-   TJ"   f (t)  e       dt    ,0).,=  211 


C   is  shown  in  Figure  3.2,  where  the   sin  x   modulation 

x 
whose  zero  crossings  with  to-  axis  are  determined  by  the 

pulsewidth. 

Example:   For  a  radar  with   PRF  of  8  00  HZ,  a  PW  of  1.3usec 
and  wavelenght  A  of  O.i  meter. 

The  first  zero  crossing  of  the  envelope  will  occur  at 
769  KHZ.   For  a  target  radial  speed  (v  )  of  900  knots, 
using  a  formula  from  ref  {6};  Doppler  frequency,  f  =  2v 


-      2  •  900  •  1852/3600      n  oc  „,7,7 

f  ,  =   rr—T =   9.2  6  KhZ 

d  O.J. 


Thus  it  is  obvious  that  the  peaks  of  the  repetitive 
"clutter  spectrum"  are  essentially  of  equal  amplitude  in  the 
frequency  range  of  interest  in  MTI-radar  (Figure  3.3). 

Figure  3.3  is  a  small  portion  of  Figure  3.2  (near  the  origin) 
and  one  moving  target  is  introduced  in  this  ideal  frequency 
spectrum  which  consists  of  impulses. 
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FIGURE  3.2. —  envelope  of  many  impulses  at 


f  =  n  •  PRF 


Figure  3.4  shows  the  "real"  clutterspectrum,  i.e.  the 

impulses  have  been  spread  out  in  frequency.   There  are 

several  reasons  for  the  nonideal  shape  shown:  The  energy 

in  the  returned  echoes  is  spread  out  in  frequence  due  to 

land  target  movement  (vegetation  swaying  in  the  wind),  sea 

waves  moving,  platform  noise  (own  irregular  movement), 

and  effects  of  a  finite  number  of  radar  returns  which  stems 

from  the  fact  that  a  radar  antenna  normally  rotates. 

Thus  we  do  not  have  an  infinite  sequence  of  target  pulses. 

Rather,  the  number  of  "hits"  on  a  particular  target  may  be, 

typically,  from  10  to  5  0  per  antenna  rotation  period. 

Figure  3.5   shows  the  bursts  of  target  echoes  every  time  the 

antenna  points  at  the  target;  the  time  between  each  burst 

is  called  the  antenna  rotation  period. 

In  Figures  3.6  and  3.7  the  Fourier  transform  of  a  finite 

number  of  6,  respectively  10  hits  are  shown.  The  mathematics 

is  omitted  here. 

The  fewer  the  number  of  hits,  the  more  energy  leaks  into 
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2PRF 


FIGURE  3.4.   A  "Real"  Clutter  Spectrum 


frequencies  where  it  may  give  raise  to  false  alarm. 
Thus  the  finite  number  of  hits  not  only  leads  to  spreading 
of  the  clutter  spectrum  around  multiples  of  PRF,  but  also 
reduces  the  maximum  attainable  cancellation  ratio.  For 
proper  MTI  operation,  at  least  10  hits  per  scan  is  needed. 
The  Fourier  transform  of  such  a  truncated  series  of  target 
returns  in  which  there  is  a  moving  target  among  heavy  clutter 
return,  is  shown  in  Figure  3.8.  Note  that  the  component  is 
also  spread  out  from  the  ideal  impulses,  and  the  energy  in 
the  sidebands. 

b.   Time  Response  Approach  to  MTI-Radar  Analysis 

In  this  analysis  we  do  not  concern  ourselves  with 
the  trailing  edge  of  a  typical  echopulse,  as  our  filter  has 
no  noticeable  disadvantageous  effect  in  this  area. 
But  since  our  filter  is  a  recursive  filter,  the  resirculation 
of  signals  lead  to  undesirable  effects  which  must  be 
investigated. 
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Video 


-1 


FIGURE  3.9.  Delay  Line  Canceller 


Figure  3.9  shows  the  blockdiagram  of  a  conventional  MTI 
delay-line  canceller  {C}.   Cancellation  is  obtained  only  when 
two  succeeding  pulses  are  equal  in  amplitude.  The  first  and 
last  pulse  in  a  pulse  train  of  a  finite  number  of  pulses  will 
be  passed  by  this  filter. 

In  the  recursive  CCD-filter  of  Figure  2.35  the  first- 
pulse  of  such  a  fixed  amplitude  pulse  train  will  be  passed, 
the  second  pulse  will  not  be  fully  cancelled.  After  the 
pulse  train  input  is  terminated,  it  will  take  a  finite  time, 
depending  upon  the  values  of  the  feedback  coefficients  b. 
and  b?,  before  the  filter  output  approaches  zero. 
This  time  is  measured  in  number  of  PRI's,  and  the  effects  in 
PPI  presentation  is  that  the  target  will  spread  out  in 
azimuth  in  the  direction  of  the  antenna  rotation. 

Figures  3.10  and  3.11  illustrates  typical  unit  impulse 
responses  of  a  first  order  recursive  filter  with  feedback 
coefficients  equal  to  0.3  and  0.9.  respectively. 
Example:  A  typical  antenna  rotation  rate  is  2  0  revolutions 
per  minute,  which  means  3  seconds/scan  or  8.3  3  msec/degree. 
For  the  case  where  b,  -  0.9  this  means  that  the  target 
would  be  extended  29/8.33  =  3.5  degrees  beyond  the  actual. 
position.  This  is  considered  to  be  a  serious  limitation, 
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ljLh(t) 


FIGURE  3.10.  Time  Response  of  First  Order  Recursive 
Filter >  b1  =  0.3 


FIGURE  3.11.  Time  Response  of  First  Order  Recursive 
Filter,  b1  =  0.9 
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and  consequently  we  have  to  keep  this  aspect  in  mind  when 
designing  filters  of  this  type. 

2 .   Applying  the  Comb  Filter  Technique. 


The  purpose  of  MTI  -radar  is  to  detect  targets 
moving  at  or  above  a  certain  minimum  speed,  and  suppress  all 
echoes  from  stationary  targets  (clutter). 
To  do  this,  one  exploits  the  well  known  doppler  effect, 
namely  that  the  frequency  of  a  returned  echo  from  such  a 


target  (moving),  is  different  from  the  transmitted  frequency. 

Using  detection  techniques  in  the  receiver  {6}  ,  {7}, only 

doppler  frequencies  are  amplified  and  become      the  MTI- 

video  which  may  be  transmitted  to  a  PPI  for  viewing. 

Thus  we  need  a  comb  filter  (Figure  3.12),  which  can  reject 

the  clutter,  and  pass  the  doppler  shifted  signals. 

One  disadvantage  is  that  when  the  target's  radial  velocity 

corresponds  to  some  multiple  of  the  PRF ,  it  will  be 

suppressed  also  (the  target  is  then  said  to  have  a  radial 

blind  speed).  To  overcome  this  one  must  vary  the  PRI  from 

pulse  to  pulse,  by  PRF-jitter  (random  variation  of  PRI), 

or  PRF  stagger  (sequentially  stepping  the  PRI  between  two 

values  only) . 

A  proposed  method  for  building  a  CCD  -  filter  to  facilitate 

such  a  system  is  outlined  in  chapter  IV. 


»> 


FIGURE  3.12.   Comb  Filter  and  Clutter 
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3  .   Bandwidth  Requirements. 

The  pulsewidth  determines  the  bandwidth  (BW) 
necessary  to  reproduce  the  signal. 


bw    2p 


The  employed  delay  lines  have  96  stages  each.  The  necessary 
clock  frequency,  f  ,  is  then: 


f   =  PRF-  9  6 
c 

For  PRF  =800  HZ,   f   =   76.8  KHZ,  which  is  the  effective 
sample  frequency  (f0). 

Using  Nyquist's  Sampling  Theorem,   f   =  2BW. 

The  maximum  bandwidth  that  can  be  faithfully  reproduced  is 
then  3  8.4  KHZ.    This  corresponds  to  PW  =  13  usee,  which  is 
quite  wide  and  unrealistic.   Thus  delay  lines  with  considerable 
larger  number  of  stages  are  needed. 

B.    EXPERIMENT 

1 .    Simulation  of  Radar  Return  Signal 

Since  the  size  of  the  delay  line  was  limited  to  9  6 
stages  it  was  decided  to  simulate  the  radar  signal  for  which 
the  delay  line  was  able  to  respond.   The  objective  was  to 
provide  radar  return  signals  simulating  fixed  and  moving 
targets  with  variable  PRF  and  number  of  hits  per  scan. 
The  simulation  of  PRF  was  implemented  using  a  Beckman  90  30 
function  generator  controlled  by  a  sinewave  generator 
(Wavetek  Model  142). 

The  function  generator  produces  a  square  pulse  whose  lenght 
and  amplitude  can  be  varied  for  every  period  of  the  sinewave. 
Thus  the  PRF  may  be  controlled  by  changing  the  Wavetek 
frequency  (Figure  3.13). 
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4*    Wavetek 


Function  Generator 


FIGURE  3.13.  Time  Relationship  Between  Output  from 

Function  Generator  and  Sinewave  Generator 


The  pulse  train  thus  obtained  is  the  simulation  of  clutter, 
or  fixed  targets,  since  each  pulse  is  equal  in  amplitude. 
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FIGURE  3.14.  Block  Diagram  of  Simulation  Configuration 
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Figure  3.14  illustrates  the  simulation  setup.   A  National 
Semiconductor  MM  4  616A  switch  (Figure  3.15)  was  used  to 
perform  the  switching  functions.  The  specifications  of  this 
arrangement  is  that  each  pulse  through  this  switch   had  a 
different  amplitude,  simulating  the  moving  target  signal. 
On  the  oscilloscope,  this  signal  should  exhibit  the  "butter- 
fly effect"  {6}.  If  needed,  both  stationary  and  moving  target 
signals  could  be  fed  into  the  filter  simultaneously. 
The  second  switch  was  operated  by  another  Wavetek  M142,  which 
was  set  such  that  the  positive  part  of  the  squarewave  was 
much  shorter  than  the  negative  part.  The  positive  part 
would  turn  the  switch  on,  letting  a  finite  number  of  pulses 
from  switch  one  through.   This  was  the  simulation  of  the 
antenna  sweep.   However,  the  present  simulation  of  this 
effect  is  limited  because  it  has  not  been  possible  to 
simulate  the  long  time  between  two  series  of  radar  hits. 

2.   Results. 


a.  Stationary  Targets  with  PRF  =  800  HZ 


FIGURE  3.17.  Cancelling  of  "Stationary"  Target, 

Horizontal:  10  psec/cm.  Vert:  .5V/cm 
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TABLE  IV 
Typical  Radar  Characteristics 


Band : 

j-band 

S-band 

C-band 

X-band   K 

-band      ' 

Frequency 
(GHZ) 

L  -  2 

2-4 

4-8 

8-12.5    1 

2.5-18 

Typical 
PRF 

200- 

500 

300- 
600 

500- 
1000 

800-      2 
3000 

000- 
5000 

PRI=  1 
PRF 

2-5 
msec 

1.67- 

3  msec 

1-2 

msec 

333- 
1250  usee 

200- 
500  usee 

Maximum 
Theoretical 
Range  n.m. 

L62- 
405 

135- 
270 

81- 
162 

27- 
101 

16.2- 
40.5 

Typical  PW 

3-5 
usee 

1-3 
usee 

0.5- 
4  usee 

.05- 
1.5  usee 

.01- 
0.5  usee 
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FIGURE  3.18. 


Cancelling  of  Stationary  Target. 
PRF=  800HZ.  Horizontal:  10  usec/cm. 


FIGURE  3.19.  Cancelling  of  Narrow,  Stationary  Target 
PRF  =  800HZ.  Horizontal:  10  usec/cm. 
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FIGURE  3.20.  Moving  Target  Signal.  PRF=  800  HZ 
Horizontal:  10  ysec/cm 


Figure  3.20  shows  a  "moving  target"  signal.  The  doppler  freq- 
uency was    at  the  center  of  first  comb  tooth  of  the  filter 
It  can  be  seen  that  the  combfilter  did  not  cancel  the  signal 
The  insertion  loss  was  quite  small. 


FIGURE  3.21.  Moving  and  Stationary  Targets.  PRF=  800  HZ 
Horizontal:  10  psec/cm 
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c.  Both  Stationary  and  Moving  Signals,  PRF  =  800  HZ 

Figure  3.21  shows  an  input  consisting  of  both  moving 
and  stationary  target  signals.  It  can  be  seen  that  only  the 
moving  signal  is  passed. 

d.  Fixed  and  Moving  Targets.  PRF  =  10  KHZ 

To  try  to  eliminate  some  of  the  glitches,  the  effect- 
ive clock  frequency  was  increased  to  28  8  KHZ,  so  that  the 
nulls  of  the  comb  filter  appeared  at  multiples  of  10  KHZ. 
The  set  of  coefficients  for  the  following  results  were 
(potentiometer  settings): 

a   =  0.6885,   an  =  0.975,   a0  =  0.87,   b,  =  0.5475,  b0=  0.491 
o  l  i  1  l 

For  all  the  pictures  the  vertical  scale  is  0.5  V/cm. 


FIGURE  3.22.  Stationary  Target  Signal.  PRF  =  10  KHZ. 
Horizontal  scale:  50  ysec/cm. 


Figure  3.22  corresponds  to  Figure  3.18,  the  difference 
being  that  now  the  PRF  was  10  KHZ  instead  of  800  HZ.  Note 
that  the  glitches   were  much  smaller  in  amplitude  at  higher 
clock  frequency. 
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FIGURE  3.23.  Moving  Target  Signal  only.  PRF  =  10  KHZ 
Horizontal  scale:  50  usec/cm. 


e.   Moving  Target.   PRF  =10  KHZ. 

Figure  3.  2  3  shows  a  moving  target  only,  with 
doppler  frequency,  f,  =  5  KHZ.  Again  it  is  noted  that  the 
moving  target  signal  passed  through  the  comb  filter  with  low 
insertion  loss. 


V 


FIGURE  3.24.  Moving  Target,  fd  =  8.5  KHZ.  PRF+  10  KHZ 
Horizontal:  50  ysec/cm/ 
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f .   Moving  Target  with  Doppler  Frequency  Close  to  the 

First  Null 

Figure  3.2  4  shows  that  the  output  amplitude  indeed 
was  attenuated  when  f ,  was  close  to  the  blind  speed. 
In  this  case  the  first  blind  speed  was  the  radial  velocity 
of  the  target  corresponding  to  f ,  =  10  KHZ. 


FIGURE  3.25.  Moving  Target  and  Stationary  Target, 
fd  =  5  KHZ,  PRF  =10  KHZ. 
Horizontal:  50  usec/cm. 

g.   Stationary  and  Moving  Targets,  PRF  =10  KHZ 

Figure  3.2  5  shows  that  when  the  input  consisted  of 
both  stationary  and  moving  targets,  only  the  moving  target 
signal  was  passed.  In  this  case  f  ,  =  5  KHZ  which  was  midband, 

h.   Stationary  and  Moving  Targets,  Doppler  Frequency  at 
First  Null 

Figure  3.26  depicts  the  same  situation  as  in  Figure 
3.25,  only  now  f  ,=  10  KHZ  which  was  the  first  blind  speed. 
Both  targets  were  suppressed,  and  the  glitches  are  seen. 
The  output  was  almost  identical  to  the  output  in  the  station- 
ary target  case  (Figure  3.22) 
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FIGURE  3.26.  Moving  and  Stationary  Targets  Superimposed 
fd  =  10  KHZ,  PRF  =  10  KHZ.  Horizontal:  50 
usec/cm. 


i.   The  Effect  of  a  Finite  Number  of  Stationary  Target 
Signals. 


FIGURE  3.27.  Truncated  Series  of  Fixed  Targets  (15  hits) 
Horizontal:  0.5  msec/cm. 

Figure  3.27  shows  the  input  and  output  of  the  comb  filter  in 
the  case  where  Switch  2  (Figure  3.14  )  was  operated  by  a  lew 
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frequency  asymmetrical  squarewave,  to  simulate  finite  number 
of  "hits"  on  the  target  during  one  period  of  the  antenna  scan 
The  input  to  the  switch  was  an  infinite  train  of  fixed 
amplitude  pulses.  The  oscilloscope  was  synchronized  to  the 
switch  controlling  waveform.   Note  that  the  first  and  last 
pulses  in  each  truncated  series  of  radar  hits  were  passed  by 
the  comb  filter.  This  was  in  agreement  with  the  time  analysis 
The  rest  of  the  pulses  were  cancelled,  except  for  the 
glitches,  which  appeared  quite  pronounced  in  this  picture. 
Additional  study  is  necessary  to  determine  the  cause  of  the 
glitches  and  their  elimination. 


FIGURE  3.28.  One  Pulse  from  Figure  3.27. 

Horizontal  scale:  10  usec/cm. 

Figure  3.2  8  shows  one  of  the  individual  pulses  in  the 
truncated  trains  of  radar  hits  from  Figure  3.27.  The 
oscilloscope  was  now  synchronized  to  the  function  generator 
(Figure  3.14).   It  is  noted  that  the  cancellation  was  poor, 
verifying  the  analysis  of  section  III. A. 1. a,  which  predicted 
that  energy  would  leak  into  sidebands  when  a  finite  number 
of  hits  were  being  processed,  causing  false  alarm. 
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FIGURE  3.29.  Truncated  Series  of  Fixed  Targets  (3  Hits). 
Horizontal:  20  usec/cm. 

Figure  3.29  shows  a  truncated  series  of  only  3  hits.  The 
oscilloscope  was  now  synchronized  to  the  squarewave  generator 
controlling  Switch  2,  and  the  three  hits  cannot  be  seen 
clearly  on  this  picture  for  that  reason.   Note  that  the  first 
two  pulses  were  passed  and  the  last  pulse  was  also  passed  as 
well  as  several  recirculated  pulses,  causing  a  smearing  of 
the  trailing  edge.  This  phenomenon  was  desribed  in  section 
III.A.l.b.,  and  is  caused  by  recirculation  of  the  signals. 
The  signal  will  die  out  gradually  depending  upon  the  size  of 
the  feedback  coefficient. 

Figure  3.30  shows  one  of  the  individual  pulses  shown  in 
Figure  3.29.   Although  this  pulse  was  similar  to  Figure  3.28, 
it  should  be  quite  different.   More  pronounced  false  alarm 
was  expected  in  Figure  3.30  than  in  Figure  3.28.  The  reason 
for  the  poor  result  here  was  that  the  individual  truncated 
series  were  too  close,  so  that  several  consecutive  series 
were  superimposed  on  the  scope.   This  was  also  the  reason 
for  the  base  line  to  appear  even  under  the  pulse  in  both 
Figures  3.28  and  3.30. 
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FIGURE  3.30.  One  pulse  from  Figure  3.29. 
Horizontal:  10  ysec/cm. 

The  relative  brightness  of  the  pulse  itself  and  the  base-line 
under  the  pulse  in  Figure  3.2  8  indicated  that  there  were  more 
oscilloscope  sweeps  for  pulses  than  for  base  line. 
But  in  Figure  3.30  the  situation  was  reversed.  This  was 
reasonable,  since  here  most  of  the  time-axis  was  a  base-line. 
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A.  CONCLUSIONS  OF  EXPERIMENTS 

1 .   Filter  Experiments 

It  has  been  shown  that  the  CTD-device  can  be  used 
for  implementing  second  order  comb  filters. 

2 .  MTI-Radar  Simulation  Experiment. 

Most  of  the  simulation  objectives  were  met  except  for 
the  simulation  of  antenna  rotation.  Due  to  the  limitation  of 
instrumentation,  the  time  between  two  series  of  radar  hits 
is  too  short.  A  more  realistic  case  would  be  to  have  a  series 
of  radar  hits  appear  every  two  seconds. 

A  major  unresolved  problem  encountered  in  this  study  was  the 
observed  glitches.  It  must  be  solved  before  the  sampled 
analog  recursive  CTD  comb  filter  may  be  successfully  employed 
in  real  radar  systems. 

Characteristics  of  the  present  CTD-devices  do  not  yet  allow 
for  enough  bandwidth.  More  delay  stages  are  needed. 
Until  then  only  simulated  results  may  be  obtained,  with 
rather  restrictive  ranges  of  pulsewidths  and  PRF ' s . 
However,  deviations  between  theory  and  measured  results  have 
been  noticed.  It  was  found  that  both  the  delay  line  and  its 
peripheral  electronics  do  not  have  adequate  bandwidths . 
Consequently,  the  coefficients  used  in  experimental  studies 
vary  with  frequency.  It  is  believed  that  this  is  one  of  the 
major  causes  of  deviation  of    experiment  from  theory. 

B.  RECOMMENDATIONS 

1 .   Research  Topics  Directly  Related  to  Present  Thesis. 

In  the  filter  study  of  chapter  II,  more  work  is 
needed  to  achieve  implementation  of  frequency  independent 
coefficients.  Perhaps  voltage  controlled  operational 
amplifiers  should  be  used  instead  of  potentiometers,  and 
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also  a  frequency  independent  delay-line  and  its  peripheral 
circuitry. 

In  the  MTI-  simulation  study,  the  problem  is  the  glitches. 
To  find  the  causes,  and  to  eliminate  them  is  important. 
Another  area  is  to  build  a  better  simulation  setup  for 
examining  the  effect  of  finite  number  of  returns. 

2.   Proposed  Application  to  Staggered/ Jittered  PRF  MTI- 
Radar . 

Not  too  many  military  MTI-radars  today  will  be 
built  that  do  not  include  either  jitter  PRF  or  stagger  PRF 
capabilities.   One  of  the  major  attractions  of  CTD-delay 
lines  is  the  ease  and  speed  with  which  the  clock-frequency, 
and  thereby  the  delay  time,  may  be  changed.   Consequently, 
the  null  frequencies  of  the  comb  filter  can  be  easily  changed 
by  changing  the  clock  frequency. 

But  there  are  major  problems  in  applying  the  CTD-filter  to 
such  cases.  This  is  because  in  jitter,  or  staggered  PRF-radar, 
the  PRI  changes  from  pulse  to  pulse,  and  thus  the  clock 
frequency  of  the  delay  line  must  be  changed  for  every  PRI. 
It  is  seen  that  information  which  has  already  been  fed  into 
the  delay-line  (memory  device) ,  will  be  clocked  out  at  a 
different  clock  frequency  and  then  added  in  the  output  summer 
(Figure  2.3)  with  a  signa   employing  a  third  PRI.  The  result 
would  be  meaningless. 

Figure  U .  1  is  a  proposal  to  employ  clocked  CTD-delay  lines 
in  a  jitter  PRF  tracking  filter  which  could  overcome  this 
problem.  It  may  solve  the  undesirable  effects  of  long 
settling  time  associated  with  high  feedback  coefficients, 
as  stated  in  section  III.A.l.b. 

A  short  explanation  using  difference  equations,  of  this 
circuit  follows: 
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1st  pulse  arrives: 


Y(0)  =  ao  X(0) 


Upper  part  of 
filter  effective: 

Video  into  lower 
circuit,  still 
tapping  upper  part 


Y(l)  =   aoX(l)  +  (b1ao+a1)x(0) 


Y(2)=aQX(2)+(b1a  +a1)X(l)  + 


bl  aox(0)  +  a1b1X(0) 


SI  opens.  Video  into 
upper  part: 


Y(3)=aQX(3)  +  (b1aQ+  a1)X(2) 


SI  closes,  S2  opens, 
Video  into  lower  part 


Y(4)=a  X(4)  +  (bna  +  an )  X(3) 
o         1  o    1 


SI  opens,  S2  closes 

Video  into  upper  part:       Y(5)  =  aQX(5)  +  (b,a  +a-L)X(U) 


etc 


If  the  video  train  input  is  a  train  of  fixed  amplitude 
pulses,  then: 


Y(n)  +  X(n){(a  (1+b. ))+an  } 
oil 


There  will  be  perfect  cancellation  if 


a  (1  +  b,  )  +  an  = 
o      11 


Circuit  operation: 
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V.   is  a  sawtooth  waveform.  A  sample  and  hold  circuit  is  used 
in  r 

to  establish  a  DC-voltage  level  which  determines  the 
frequency  output  of  the  voltage  controlled  frequency  genera- 
tors. Starting  with  frequency  generator  #  1,  assuming  that 
the  video  is  being  fed  into  the  upper  part  of  the  filter,  the 
binary  counter  will  count  9  6  periods  and  then  generate  a 
trigger  pulse,  which  controls  the  radar  transmitter  and  all 
radar  timing  circuits  and  the  various  switches  in  the  filter 
itself.  The  functions  of  this  trigger  pulse  in  the  filter  are 
to: 

-  Reset  the  lower  counter  to  zero. 

-  Sample  V .   to  establish  a  new  f  0  . 

r     in  c2 

-  Switch  the  video  to  the  lower  part  of  the  filter. 

-  Switch  off  S2,  the  output  of  the  lower  delay  line. 

-  Switch  on  SI,  the  output  of  the  upper  delay  line. 

Next,  after  another  PRI  and  f  «  counted  to  96,  a  trigger  pulse 

appears  (the  counters  will  only  produce  a  trigger  pulse  on  the 

count  of  9  6  after  each  reset),  f  -  is  then  set  to  a  new  value, 

cl 

counter  #  1  is  reset,  the  video  is  switched  to  the  upper  part 
of  the  filter,  the  switch  SI  is  opened  and  S2  is  closed. 
Thus  the  direct  signal  will  always  be  summed  in  one  part  of 
the  filter  with  the  delayed  signal  stored  in  the  other  part. 
This  scheme  may  be  extended  to  higher  order  filters. 
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APPENDIX  A 

FILTER  DESIGN  IN  THE  Z  -  DOMAIN 

It  was  desired  to  design  a  filter  with  100  HZ  across 
the  null  at  an  attenuation  of  at  least  2  0  db  down. 
From  a  pole-zero  plot  in  the  z-domain  Figure  A-l),  the 
magnitude  of  the  transfer  function  can  be  evaluated  for  every 
point  on  the  unit  circle,  geometrically. 
By  changing  the  pole  and  zero  locations  within  the  unit 
circle  in  a  trial  and  error  manner,  and  comparing  the 
magnitude  of  the  transfer  function  with  the  desired  shape  of 
the  frequency  response  curve,  it  was  possible  to  come  fairly 
close  after  a  rather  short  time: 


Z-plane 


FIGURE  A-l.  Pole-Zero  Plot  in  Z-Domain. 


Results:   Zeros:  Z12  =  0.8  ±  jO.l 
Poles:  P12  =  -0.1  ±  jO.5 


This  gives  the  transfer  function: 


H(z) 


0.556  -  z"1  +  0.U556  z"2 
1  +  0.14  z"1  +  0.2074  z"2 
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A  computer  simulation  of  a  filter  using  these  coff icients 
is  shown  in  Figure  A-2.  The  desired  width  of  the  stopband , 
and  attenuation  was  obtained.  The  filter  was  also  built, 
and  a  few  data  points  were  plotted,  but  the  implementation 
of  the  coefficients  were  not  correct. 
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